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Abstract 

We compute the fully non-linear Cosmic Microwave Background (CMB) anisotropies on scales 
larger than the horizon at last-scattering in terms of only the curvature perturbation, providing a 
generalization of the linear Sachs- Wolfe effect at any order in perturbation theory. We show how 
to compute the n-point connected correlation functions of the large-scale CMB anisotropies for 
generic primordial seeds provided by standard slow-roll inflation as well as the curvaton and other 
scenarios for the generation of cosmological perturbations. As an application of our formalism, we 
compute the three- and four-point connected correlation functions whose detection in future CMB 
experiments might be used to assess the level of primordial non-Gaussianity, giving the theoretical 
predictions for the parameters of quadratic and cubic non-linearities /nl and c/nl- 

PACS numbers: 98.80.cq 



* Electronic address: 
^Electronic address: 

* Electronic address: 



nbartolo@ictp.trieste.it 

sabino.matarrese@pd.infn.it 

antonio.riotto@pd.infn.it 



I. INTRODUCTION 

Cosmological inflation [l] is the dominant paradigm to understand the initial conditions 
for CMB anisotropies and structure formation. In the inflationary picture, the primordial 
cosmological perturbations are created from quantum fluctuations "redshifted" out of the 
horizon during an early period of superluminal expansion of the universe, where they remain 
"frozen" . They are observable as temperature anisotropies in the CMB at the last scattering 
surface. They were first detected by the Cosmic Background Explorer (COBE) satellite 
^| . The last and most impressive confirmation of the inflationary paradigm has been recently 
provided by the data of the Wilkinson Microwave Anisotropy Probe (WMAP) mission |5[. 
Since the observed cosmological perturbations are of the order of 10 -5 , one might think that 
first-order perturbation theory will be adequate for all comparison with observations. That 
may not be the case however, as the Planck satellite p and its successors may be sensitive to 
the non-Gaussianity of the cosmological perturbations at the level of second- or higher-order 
perturbation theory Statistics like the bispectrum and the trispectrum of the CMB can 
be used to assess the level of primordial non-Gaussianity on various cosmological scales and 
;o discriminate it from the one induced by secondary anisotropies and systematic effects 



7, 



1CJ. Therefore, it is of fundamental importance to provide accurate theoretical 
predictions for the statistics of the large-angle CMB anisotropies as left imprinted by the 
primordial seeds originated during or immediately after inflation. Steps towards this goal 
have been taken in Refs. H, 12,Q,Q| at the level of second order perturbation theory. 

In this paper we derive an expression for the anisotropies of the CMB on scales larger than 
the horizon at last scattering which is valid at any order in perturbation theory, providing a 
fully non- linear generalization to the Sachs- Wolfe effect at first- [15] and second-order |l3| . 
In particular, for the standard single-field models of inflation, we provide the exact non- 
perturbative expression for the bispectrum and the trispecturm in the so-called "squeezed" 
limit in which some of the wavenumbers are much smaller than the others. Furthermore, we 
compute the generic expressions for the non-linearity parameters /nl and #nl characterizing 
respectively the quadratic and cubic non-linearity in the large-angle CMB anisotropies. 

The paper is organized as follows. In Section II we provide the non-linear generalization 
of the Sachs- Wolfe effect which is expressed in terms of the comoving curvature perturbation 
from inflation in section III. In Section IV we show how to compute the connected n-point 
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correlation functions of the CMB anisotropies, leaving the details to the Appendix. Finally, 
we draw our conclusions in Section V. 



II. THE NON-LINEAR SACHS-WOLFE EFFECT 

Our starting point is the Arnowitt-Deser-Misner (ADM) formalism which is particularly 
useful to deal with the non-linear evolution of cosmological perturbations. The line element 
is 

ds 2 = —N 2 dt 2 + Ni dt dx i + dx i dx j , (1) 

where the three-metric 7^, the lapse N and the shift Ni functions describe the evolution 
of timelike hypersurfaces. In the ADM formalism the equations simplify considerably if we 
set N l = 0. Moreover we are interested only in scalar perturbations in a flat Universe and 
therefore we find it convenient to recast the metric as 

d s 2 = - e 2 * dt 2 + a 2 (£)e" 2 % dx i dx j , (2) 

where a(t) is the scale factor describing the evolution of the homogeneous and isotropic Uni- 
verse and we have introduced two gravitational potentials <3? and The expression (J2J) holds 
at any order in perturbation theory. To make contact with the usual perturbative approach, 
one may expand the gravitational potentials at first- and second-order, e.g, $ = $1 + $2/2. 
From Eq. (J2J) one recovers at linear order the well-known longitudinal gauge, N 2 = (1 + 2$i) 
and 7ij = a 2 (l — 2^1)6^. At second-order, one finds $2 = 02 — 20 2 and ^2 = "02 + 2 , 2 where 
0i, -01 an d 02, "02 (with 0i = $1 and ipi = ^1) are the first and second-order gravitational 
potentials in the longitudinal (Poisson) gauge adopted in Refs. N 2 = (1 + 20i + 2 ) 

and jij = a 2 (l — 2^i — ■0 2 )^ij as far as scalar perturbations are concerned. In writing Eq. (j2J) 
we have neglected vector and tensor perturbation modes. For the vector perturbations the 
reason is that we are interested in long- wavelength perturbations, i.e. on scales larger than 
the horizon at last scattering, while vector modes will contain gradient terms being pro- 
duced as non-linear combination of scalar-modes and thus they will be more important on 
small scales (linear vector modes are not generated in standard mechanisms for cosmological 
perturbations, as inflation). For example the results of Ref. show clearly this for second- 
order perturbations. The tensor contribution can be neglected for two reasons. First, the 
tensor perturbations produced from inflation on large scales give a negligible contribution 



to the higher-order statistics of the Sachs- Wolfe effect being of the order of (powers of) 
the slow-roll parameters during inflation (this holds for linear tensor modes as well as for 
tensor modes generated by the non-linear evolution of scalar perturbations during inflation, 
for example see the results of Ref. [24] for second-order perturbations). Moreover, while on 
large scales the tensor modes have been proven to remain constant in time 17], when they 
approach the horizon they have a wavelike contribution which oscillates with decreasing 
amplitude. 

Since we are interested in the cosmological perturbations on large scales, that is in per- 
turbations whose wavelength is larger than the Hubble radius at last scattering, a local 
observer would see them in the form of a classical - possibly time-dependent - (nearly zero- 
momentum) homogeneous and isotropic background. Therefore, it should be possible to 
perform a change of coordinates in such a way as to absorb the super-Hubble modes and 
work with a metric of an homogeneous and isotropic Universe (plus, of course, cosmological 
perturbations on scale smaller than the horizon). We split the gravitational potential $ as 



where $£ stands for the part of the gravitational potential receiving contributions only from 
the super-Hubble modes; $ s receives contributions only from the sub-horizon modes 



where H is the Hubble rate computed with respect to the cosmic time, H = a/a, and 8(x) 
is the step function. Analogous definitions hold for the other gravitational potential 

By construction $^ and are a collection of Fourier modes whose wavelengths are larger 
than the horizon length and we may safely neglect their spatial gradients. Therefore and 

are only functions of time. This amounts to saying that we can absorb the large-scale 
perturbations in the metric (j2J) by the following redefinitions 



(3) 




(4) 



dt 



e*<dt , 



(5) 



a = ae 



(6) 
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The new metric describes a homogeneous and isotropic Universe 



ds 2 = —dt 2 + a 2 Sij dx l dx-' , 



(7) 



where for simplicity we have not included the sub-horizon modes. On super-horizon scales 
one can regard the Universe as a collection of regions of size of the Hubble radius evolving 
like unperturbed patches with metric (|ZJ) [121 • 

Let us now go back to the quantity we are interested in, namely the anisotropics of the 
CMB as measured today by an observer O. If she/he is interested in the CMB anisotropies 
at large scales, the effect of super-Hubble modes is encoded in the metric (JJJ). During their 
travel from the last scattering surface - to be considered as the emitter point £ — to the 
observer, the CMB photons suffer a redshift determined by the ratio of the emitted frequency 
Us to the observed one uJo 



where Tq and T% are the temperatures at the observer point and at the last scattering 
surface, respectively. 

What is then the temperature anisotropy measured by the observer? The expression (JSJ) 
shows that the measured large-scale anisotropies are made of two contributions: the intrinsic 
inhomogeneities in the temperature at the last scattering surface and the inhomogeneities 
in the scaling factor provided by the ratio of the frequencies of the photons at the departure 
and arrival points. Let us first consider the second contribution. As the frequency of the 
photon is the inverse of a time period, we get immediately the fully non-linear relation 



As for the temperature anisotropies coming from the intrinsic temperature fluctuation at 
the emission point, it maybe worth to recall how to obtain this quantity in the longitudinal 
gauge at first order. By expanding the photon energy density p 7 oc the intrinsic temper- 
ature anisotropies at last scattering are given by SiTg/Ts = (l/4)5ip 7 /p 7 . One relates the 
photon energy density fluctuation to the gravitational perturbation first by implementing 
the adiabaticity condition 5ip 7 /p 7 = (4/3)<5ip m /p m , where Sip m /p m is the relative fluctua- 
tion in the matter component, and then using the energy constraint of Einstein equations 
$1 = — (l/2)5ip m /p m . The result is 5iT £ /T £ = — 2$ 1£ /3. Summing this contribution to 
the anisotropies coming from the redshift factor (|UJ) expanded at first order provides the 




(8) 



(9) 
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standard (linear) Sachs- Wolfe effect b{FojTo = $if/3. Following the same steps, we may 
easily obtain its full non-linear generalization. 

Let us first relate the photon energy density p 7 to the energy density of the non-relativistic 
matter p m by using the adiabaticity conditon. Again here a bar indicates that we are 
considering quantities in the locally homogeneous Universe described by the metric 
Using the energy continuity equation on large scales dp/dt = —3H(j> + P), where H = 
dlna/dt and P is the pressure of the fluid, one can easily show that there exists a conserved 
quantity in time at any order in perturbation theory jl9| 

T = \na + \( P ,_ dp . . (10) 
3 J (p' + P'' 



The perturbation bT is a gauge-invariant quantity representing the non-linear extension 
of the curvature perturbation ( on uniform energy density hyper surf aces on superhorizon 
scales for adiabatic fluids Indeed, expanding it at first and second order one gets the 
corresponding definition (i = —ipi — <5ip/p and the quantity ( 2 introduced in Ref. |2p| . 
At first order the adiabaticity condition corresponds to set d 7 = Cim f° r the curvature 
perturbations relative to each component. At the non-linear level the adiabaticity condition 
generalizes to 

1 f d Pm = 1 f dp-, 

3J Pm 4J p 

or 



SJ p m A J p 7 



lnp m = In pf. (12) 

To make contact with the standard second-order result, we may expand in Eq. ()12)) the 
photon energy density perturbations as <5p 7 /p 7 = 5ip 7 /p 7 + |5 2 p 7 /p 7 , and similarly for the 
matter component. We immediately recover the adiabaticity condition 

5 2 p~t _ 4: 5 2 p m | 4 / £ip m \ 2 ^ l3 ^ 
!_. p 7 3 p m 9 \ p m J 

given in Ref. [7|. 

Next we need to relate the photon energy density to the gravitational potentials at the 
non-linear level. The energy constraint inferred from the (0-0) component of Einstein equa- 
tions in the matter-dominated era with the "barred" metric (J7|) is 

—2 8irG N _ 

H =^^P m - (14) 
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Using Eqs. (jSJ) and (0) the Hubble parameter H reads 



-jj Ida 

a dt 



(15) 



where H = d In a/dt is the Hubble parameter in the "unbarred" metric. Eq. (fPfj) thus yields 
an expression for the energy density of the non-relativistic matter which is fully nonlinear, 
being expressed in terms of the gravitational potential $^ 



Pm ~ Pm& 



-2<S> f 



(16) 



where we have dropped ^£ which is negligible on large scales. By perturbing the expression 
f!16|) we are able to recover in a straightforward way the solutions of the (0-0) component 
of Einstein equations for a matter-dominated Universe in the large-scale limit obtained at 
second-order in perturbation theory. Indeed, recalling that $ is perturbatively related to the 
quantity _ <t> = <fii +02/2 used in Ref. ^ by <3>i = 0i and $2 = 4>: 



obtains 
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21] 



2(0i) 2 , one immediately 



$lPm , 

= -20i, 

Pm 

Y— = -0 2 + 4(0i) 2 - (17) 

^ Pm 

The expression for the intrinsic temperature of the photons at the last scattering surface 
Ts oc pi/ 4 follows from Eqs. (|T2^) and (JTHJ) 

T £ = T £ e- 2 ^l\ (18) 

Plugging Eqs. and (fTSj) into the expression (jHJ) we are finally able to provide the ex- 
pression for the CMB temperature which is fully nonlinear and takes into account both the 
gravitational redshift of the photons due to the metric perturbations at last scattering and 
the intrinsic temperature anisotropies 

To={^)T e e^\ (19) 

From Eq. (|19j) we read the non-perturbative anisotropy corresponding to the Sachs- Wolfe 
effect 

;*'/ 3 -l. (20) 



SnpTo 



To 



Eq. ()20|) is one of the main results of this paper and represents at any order in perturbation 
theory the extension of the linear Sachs- Wolfe effect. At first order one gets 

8 X T 



To 



V 

3 ' 



and at second order 



16*T. 



2J-0 



2 Tt 



o 



6 18 V ; 



(21) 



(22) 



which exactly reproduces the generalization of the Sachs- Wolfe effect at second-order in the 



perturbations found in Ref. 



7. 



2l( (where $i = fa and $ 2 = <p2 - 2(</>i) 2 ). 



III. RELATING THE CMB ANISOTROPIES TO THE INFLATIONARY CO- 
MOVING CURVATURE PERTURBATION 

In this section we relate the gravitational potentials <J>£ to ^ £ to the curvature perturbation 
( = bT at any order in perturbation theory (for notational simplicity we drop the subscrip 
from now on). This will allow to express the non-linear temperature fluctuations in 
terms of the initial conditions provided by inflation. We use the evolution equation in the 
ADM formalism 

- N\ + l -N\5\ + N = N SvrG S\ , (23) 

where a vertical bar denotes a covariant derivative with respect to 7^ and we have used that 
fact that the traceless part of the extrinsic curvature Kij = K{j — Kjij/3 vanishes in our 
metric (jJJ). The extrinsic curvature is defined as = —j'ij/(2N) and K = j* 3 Kij. Here 
= Y j ®Rjk =® R\ R where ®Rjk is the Ricci tensor of constant time 
hypersurfaces associated with the metric 7^. Analogous definitions hold for S % k constructed 
from the matter stress-energy three-tensor Sij = T^. 

From Eq. (|23|) we want to obtain a constraint between the gravitational potentials $ and 
^ keeping track of non-local terms. Thus we cannot neglect the gradient terms appearing 
in this equation. Since = we find 

e e 



2* 



2* 

{3) i? = 2—- (2V 2 ^ - y^' 1 ) . (24) 

(X I Xi J 



Hence Eq. (J23)l reads 

+ (V fe + V 2 ^5\, + - * m #' m ^fc) - jj(2V 2 ^ - tf,!*' 1 )^ = 87rGa 2 (t)e- 2 *^ fc , 

(25) 

where the indices of the partial derivatives are raised by 5 U . Notice that at first order this 
equation gives the usual constraint $1 = As far as the matter content is concerned we 
can consider the perfect fluid energy momentum tensor T^ v = (j> + p\ u^u v + P g^ v and we 
find 



S\ = Y% k = a" 2 (t)e 2 * (p + P) u'uk + P5\ , 
S i k = a - 2 (t)e 2vE ' (p + T) (u'uk - iju^JZ) , 



(26) 



where u % = S^Uj. We need an expression for the spatial velocities. Thus we use the 



momentum constraint of the ADM equations which reads 

8%G Ji 



22 



Ja-„, 



(27) 



where Ji = NT^ = e*T°j = — e~*Toi is the momentum density and K = —3e (H(t) — 
\&). One can use the normalization g^u^Uy = —1 to express uq, finding uq = —(1 + 
a~ 2 (t)e 2 * V-Uj) 1 / 2 . Let us just consider scalar velocities u % = d^u (on sufficiently large scales 
vector modes may be safely neglected). From Eq. (|27jl we obtain 

(p + P) Ui (l + a~ 2 (t)e 2 V^) 1 / 2 = [^(H(t) - ¥)] ! . 

In fact we are interested in the case of non-relativistic matter P m = and p m = p m e~ 
As Eq. (f25j) contains at least two gradients (also S k contains at least two spatial gradients), 
using the gradient expansion we may restrict ourselves to the solution 

^ - - e*H(t) 



(28) 



-2$ 



Ui 



AuGp 

where we have neglected We thus find 



4:7rG p r 



(29) 



2* 



Q 1 

^ k 



6nG a 2 



(30) 



Inserting Eq. (j30j) into Eq (|25|) and applying the operator did k we find 

V 4 (W _ $) = (V* fe ) + ^V 2 

+ 3d t d k fc ) - V 2 . (31) 

Notice that Eq. (}3*T]) is the non-linear generalization of the constraint between the grav- 
itational potentials $ and \I> in the longitudinal (Poisson) gauge and is valid at any or- 
der in perturbation theory. At linear order one recovers the well-known result $1 = ^i, 
while at second order the relation first found in Refs. follows 1p2 ~ <p2 = + 

10V-%d k (^i,k) ~ f V- 2 (^iVi,i) (where $ 2 = 2 - 20? and ^2 = ^2 + ^\ and fa = <f> u 
fa = 

In the following we find it convenient to write the relation (}3*Tj) as ^ = $ + /C[$, where 
the kernel /C is obtained by acting through the operator V -4 onto the r.h.s. of Eq. ()31)1 and 
takes into account non-local terms. Going to momentum space, one easily realizes that in 
the specific "squeezed" limit of Ref . 2J| , where one of the wavenumbers is much smaller than 
the other two, e.g. k\ <C ^,3, the kernel fC — > 0. From Eq. (fTUj) the curvature perturbation 
is given by 

C = ^=-*+iln^ = -*-^, (32) 
op 3 

where we have used p m = p m e~ 2<l '. Hence we finally find — 5$/3 = ( + JC and we can recast 
the non-linear temperature anisotropies (J2(J|) as 

%^ = e-^/s - 1 , (33) 
To 

which is the starting point to evaluate the n— point correlation function for CMB temper- 
ature anisotropies. The comoving curvature perturbation ( is conserved at any order in 
perturbation theory on large scales and therefore one can fix its properties right at the end 
of inflation. Since for standard inflation the curvature perturbation can be considered as a 
Gaussian distributed quantity (deviations from non-Gaussianity are proportional to 

deviations of the spectral index from unity and are, therefore, tiny), we will adopt ( as our 
Gaussian seed in the case of standard single-field models of inflation. 
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IV. CORRELATION FUNCTIONS OF LARGE-SCALE CMB ANISOTROPIES 



As a warm-up excercise, let us first consider a simpler case and evaluate the n-point cor- 
relation function of e v< ^> where <£>(x) is a Gaussian random field. Applying the well-known 
techniques of quantum field theory, it turns out that 

( e P( x e v( x iv)^ — e hl ^xdyJ(x) (<p(x)tp(y)) J(y) ^ ^34^ 

where J(x) = E^Li <5( x — x «) corresponds in fact to the source term appearing in the path 
integral formulation. The calculation of the integral in Eq. (JMj) brings 

( e v(xi) e ^(xjv)^ _ e IE;,jM x iM x i)) _ (3^ 

Notice, for example, that for the 2-point correlation function we find the usual result 

^g^( x l) e ^( x 2)^ — e if 2 ) gM x lM x 2)) ^ ^3g^ 

which for xj = x 2 gives (e 2ip( -^) = e 2 ^ 2 \ where (y? 2 ) = (y? 2 (x)). 
For the 3-point function one finds 

where for simplicity we have used the notation ((fifj) = (^(xj)yj(x J -)). If now we expand 
the exponential in the limit in which the two point function is small we obtain 



^ (xi) e^ (x2) e v(x3) ) ~ e* {v > [l + ((^ 2 ) + cycl. 

+ ^ ((^i^s) 2 + cycl.) 

+ ((^2) (^3) + cycl.) 1 . (38) 



It is in fact a term analogous to the last contribution in Eq. ()38|) that enters in the three-point 
function of the CMB anisotropies (J2"U|) if we are interested in non-linearities up to second- 
order terms only. The complication arises from the fact that the gravitational potential $ 
appearing in Eq. (J2Uj) is not a Gaussian variable, since it can already contain quadratic (and 
higher order) terms in the Gaussian variable ( in the case of single-field models of inflation. 
These will add non-Gaussian contributions which are contained in the kernel /C. In other 
scenarios for the generation of the cosmological perturbations ( is not a Gaussian quantity 
and we will properly take into account the primordial non-Gaussian contributions. 
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Let us see how to generalize the previous procedure to this case. First of all we apply 
an iterative procedure to express the kernel /C in terms of powers of (. We will use Eqs. 
(131)) and (|32|) . For the zeroth and first-order terms of the iteration we find (the suffix 
does not refer to the order of the expansion in the perturbations, but to the order of the 
approximation given by the iteration procedure: each r-th term contains up to (r + 1) powers 

of C) 

$(°) = ~C (39) 
5 

= 4C-(T) 3 /C[C 2 ], (40) 



5" V5 

and for the next terms (n = 1,2,..) 



$(2n) = $(2n-l) + ^ [$(°) ; $(2™~2) _ $ (2n-l)j ( 4i ) 
n-2 

_l_ A^ 1 [$( ,n ) — $( m+1 ) ; (j>(2n-m-3) _ gj(2n-m-2) j ^ 

m=0 

$(2n+l) = $ (2n) _j_ ^ 1 [$(0) ) $(2n-l) _ $ (2n)j 
n-2 

_|_ A^ 1 [$( m ) — $( m+1 ) ) $(2n-m-2) _ ^(2n-m-l)j 

m=0 

+ /C 2 [($ (n_1) - $ (n) ) 2 ] , (42) 



where we have introduced the bilinear operators 



£i[(0,(01 = v- 4 (6^[(-)' 4 (-),.]-2V 2 [(-)' 4 (-),]) 
/c 2 [(-),(-)] = 



1 



v- 1 i^d t d k 



(43) 



Notice that for equal entries /Ci = 6/C/5. If the upper limit of the sums appearing in these 
expressions turns out to be negative the sum must be taken to be vanishing. 

Thus we can use Eq. (|39|) and ()4Uj) to find the next order approximations for the expression 
of $ in terms of the Gaussian curvature perturbation £. For example for $( 2 ) up to 0(C 3 ) 
contributions we have 

,3,3 



* (2) = -f c_ GD nc\+^ 



(44) 



Using this iterative procedure we express the kernel K, in Eq. (|33|) as a function of £ 
and we can determine the n-point connected correlation functions for the temperature 
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anisotropies (|3*3^) accounting for the information contained in JC[Q. In the following we 
outline the procedure and we give the results, while more details on the computation can 
be found in the Appendix. 

In order to obtain the n-point (connected) correlation functions we will borrow again some 



techniques of functional-integral analysis from quantum field theory 125 
applications of the path-integral approach in cosmology, see Refs. [28 



26 



29 



30 



271 



31 



or different 
32, 33 V In 



particular we need the so called generating functional of the correlation functions W, which 
in our case can be written as 



Z[J} = J P[C]P[C]e^ dxJ(x)(e 



-C/5-K(C)/5_ 



1) 



(45) 



Here J(x) is an external source perturbing the underlying statistics, V[C] is the probability 
density functional which in our case is a Gaussian one (see Eq. (jA.3|) of the Appendix), and 
T>[Q is a suitable measure such that the total probability is normalized to unity, / 2}[C]^[C] = 
1. The kernel JC[Q will be given by the iterative procedure. From the generating functional 
one can obtain the correlation functions by taking the functional derivatives of Z\J] with 

n 

respect to the source J evaluated at J = |26J 

-Cn/5-K(Cn)/5 



^ e -Ci/5-/C(Ci)/5 _ 1 y 
,-n 6 n Z[J] 



1)) 



5J(xi)...5J(x n ) 



Z (n) { 



Xl 



(46) 



(J=0) 



where Q = C( x i)- 

By defining the new functional W [J] = lnZ[J] one can obtain the connected correlation 
functions by the same operations 



// e -Ci/5-x:(Ci)/5 . 
,-n S n W[J] 



IV. 



,-C»/B->C(C»)/5 



l)>c 



6J(x 1 )...6J(x n ) 



(47) 



(J=0) 



It is clear that the complete knowledge of the statistical properties of the perturbations, 
i.e. the complete knowledge of the correlation functions at all orders, can be achieved if one 
knows the generating functionals. In turn from the definitions ()46|) and (}4Tj) the correlation 
functions will appear in a power series expansion of the generating functionals, as 



Z[J] 
W[J] 



oo -n ~ 

1 + E - / dx 1 ...dx n Z^(x 1 , ...,x n ) J( Xl )...J(x r , 

n=l U - J 

CO jfl . 

-r / dx 1 ...dx n W( n \x 1 , ...,x n ) J( Xl )... J(x n ) . 



(48) 
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To compute the functional derivatives of W[J] one can follow the standard procedure 
to evaluate the connected correlation functions used in field theory [25]. The computation 
makes use of a perturbative expansion around some known solution, which corresponds to 
the connected correlation functions of the free scalar field. In our case, the known solu- 
tion corresponds to the correlation functions when the kernel K, vanishes, which we have 
computed previously in Eq. (|35|). As explained in detail in the Appendix, we perform the 
expansion around /C = using as expansion parameter the r.m.s amplitude of the cosmo- 
logical perturbations themselves, that is (C 2 ) 1 ^ 2 ~ 10~ 5 . 



A. The Bispectrum 

From the general results presented in the Appendix we provide now the expression for the 
3-point connected correlation function. It suffices to expand the kernel /C up to second order. 
The total kernel can be written as a convolution in configuration space 

/C(C) = J dxidx 2 tf 2 (x - xi, x - x 2 )C(x 1 )C(x 2 ) , (49) 

where K(~x — x x , x — x 2 ) is the double inverse Fourier transform of the expression 

K 2 (k 1; k 2 ) = (a NL - 1) + I Uk, • k 3 )(k 2 • k 3 )A 4 - (l/3)(k 1 • k 2 )/k 2 } , (50) 

5 L J 

where k = |k 3 | and k 3 = — (ki + k 2 ). In Eq. (|4T)|) we have added the constant onl whose 
role is to parametrize the primordial non-Gaussianity generated during or after inflation 
in the various scenarios for the generation of the cosmological perturbations, C = Cl + 
(&nl — 1) * Cl (fro m n ow on we will remove the subscript "L"). For instance, in single field 
models of inflation cznl = 1 (plus tiny contributions proportional to the deviation from scale 
invariance), in the curvaton scenarion cinl = (3/4r) — r/2, where r ~ (p a / ' p)o represents the 
relative curvaton contribution to the total energy density at curvaton decay Q|. 
We find 

iy( 3 )( Xl ,x 2 ,x 3 ) = ((e-«^)/ 5 - l)(e-^)/ 5 - l)(e-^)/5 _ l)) connected 
- 5^y dy 1 rfy 2 A" 2 (x P2 -yi,x P2 - y 2 ) [w 2 (x pi , y x ) 



x w 2 (x P3 ,y 2 ) + -w 4 (Xp 1 ,Xp 3 ,y 1 ,y 2/ 
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(51) 



where the sum is over all permutations Pi,P2,P3 taking the values (1, 2, 3) and 

uS 2 (x 1)X2 ) = - l -({e-^ - l)C 2 ))co„ n . = -^e« 2 >/ 50 (CiC 2 ) , (52) 

^ 4 (X1,X 2 ,X 3 ,X 4 ) = 4((e- Cl/5 - l)(e~^ 5 - l)C 3 C4)conn. (53) 



25 

e (C 2 >/25 

3 x 5 4 



( e (C 1 C 2 >/25 _ ^ ((CiC4) + (CaC4)) ((CiCg) + (CaC3)) + cydic 



Despite the fact that we expanded the kernel up to second-order, the expression ()51|) 
becomes exact at any order in perturbation theory in the squezeed limit for which K. tends 
to zero and for single field models of inflation for which onl = 1. In such a case, the exact 
three-point correlation function for the temperature anisotropics on large-scales is 

H/( 3 )(x 1 ,x 2 ,x 3 ) = iy (2) (x 1 ,x 2 )^ (2) (x 1 ,X3) + iy (2) (x 1 ,X3)iy (2) (x 2 ,x 3 ) (54) 

+ iy (2) ( Xl ,x 2 )^ (2) (x 2 ,x3) + iy (2) (x 1 ,x 2 )iy (2) (x 2 ,x 3 )iy (2) (x 3 ,x 1 ) , 



where 



(2k)* 

J^joiklxi-XjDVdk), (55) 



1 f d In k 
50 



and Vc, = A(ko) 2 (k/ko) ns ~ 1 is the primordial power spectrum of the comoving curvature per- 
turbation with amplitude A and spectral index rig- The expression for the exact bispectrum 
of temperature anisotropies valid at any order in perturbation theory is 

^M^^ = (2nf ( kl + k2 + k3 ) B (kx, k 2 , k 3 ) , (56) 
where 



S(k 1; k 2 ,k 3 ) = P(h)P(k 2 ) + P(h)P(k 3 ) + P(k 2 )P(h) 

+ /^Pdq-k^Pdq-kaDPdq-kal) (57) 

~ 2P(h)P(k 2 ) + J -^P(|q|)P(| q _k 2 |)P(|q + k 2 |) , < fc 2 , fc 3 ) . 

For generic momenta configurations and for models for which onl is sizeable, the exponentials 
present in the first line of Eq. (foTj) has to be expanded at second-order to consistently match 
the order of the kernel K 
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-C/5 



5' ' 2x5- (58) 
We immediately recover the expression obtained at second-order in perturbation theory for 

the non-linearity parameter /nl defined as the coefficient of the gravitational potential $ 

expanded at second-order in terms of the linear Gaussian gravitational potential $l = — <f>i, 

$ = $l + /nl * (^l) 2 , with the convential Sachs- Wolfe effect expressed as 5T/T = — ($/3) 

m 



NL 



-(l-a NL ) + - 



(59) 



+ [3(k 1 -k 3 )(k 2 -k 3 )A 4 -(k 1 -k 2 )A 2 . 

It may be worth noticing that the coefficient 1/6 in the last expression is simply the result 
of the expansion of the exponential ()58|) expressed in terms of the gravitational potential 
X : e ~C/ 5 = e ~® L/3 ~ 1 - |$ L + which gives the contribution 3 x (1/18) = 1/6 to / NL . 



B. The Trispectrum 



Let us now follow a similar procedure to obtain the 4-point connected correlation function. 
In this case the kernel K, appearing in Eq. ()33|) must be expanded up to third order. In 
configuration space it can be written as a convolution 

K(C) = J dyL 1 dy. 2 K 2 {* - x l5 x - x 2 )((xi)((x 2 ) 

+ y<ixi(ix 2 (ix3ir3(x-Xi,x-x 2 ,x-X3)C(xi)C(x 2 )C(x3), (60) 

where K 2 is the kernel defined by Eqs. (I49|) and Eq. ()50|). while K 3 (-x — x 1; x — x 2 , x — x 3 ) 
is the triple inverse Fourier transform of the expression 

#3(k 1; k 2 , k 3 ) = (6 NL - 1) + (a NL - 1) A(kx, k 2 , k 3 ) + C(k 1; k 2 , k 3 ) , (61) 

with 



^4(ki,k 2 ,k 3 



6 



ki ■ k 4 (k 3 • k 4 + k 2 ■ k 4 ) + (k 2 ■ k 4 )(k 3 ■ k 4 ) 

k A 



1 kx ■ (k 2 + k 3 ) + k 2 • k 3 

3 



k 2 



C(k!,k 2 ,k3) 



54 (k 4 ■ k 3 ) [(ki + k 2 ) ■ k 4 ] 
25 A: 4 
1 ki-k 2 



(k 1 -(k 1 + k 2 ))(k 2 -(k 1 + k 2 )) 
|ki + k 2 | 4 



3 Ikx + k 2 | 2 



+ cycl. 



(62) 



(63) 
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where k = |k 4 | and k4 = — (ki + k2 + k 3 ). In Eq. ()63|) one has to take cyclic terms by an 
exchange of the wavenumbers ki, k2, k3. 

In order to compute K 3 (ki, k 2 , k 3 ) we have applied the iterative procedure described 
in Section IV taking into account also a possible primordial non-Gaussian contribution by 
expanding the curvature perturbation as 

C = CL + (aNL-lWL + (frNL-l)*CL- (64) 

The value of &nl will depend on the different scenarios for the generation of the cosmological 
perturbations. For example, for standard single-field models of inflation &nl = 1 (plus tiny 
contributions proportional to powers of the slow-roll parameters), while for other scenarios 
it might well be non- negligible. For simplicity from now on we will remove the subscript 
"L". 

Similarly to the bispectrum, also for the 4-point connected correlation function there 
exists a specific limit for which the kernel /C tends to zero, corresponding to take two 
wavenumbers much smaller than the other ones, e.g. k\, fc 2 *C ks, ki. For this limit and in 
the case of single-field models of inflation (a NL = 1,&nl = 1), one can compute an exact 
expression of the trispectrum by Fourier transforming the 4-point connected correlation 
function which, in this limit, reads 

iy( 4 )(x 1 ,x 2 ,x 3 ,x 4 ) = ((e-^ 5 - l)(e-^ 5 - l)(e-^ 5 - l)( e -^ 5 - l)) conn . 
= iy (2) (xi, x 2 )W (2) (x 1 , x 3 )W (2) (x 1) x 4 )iy (2) (x 2 , x 3 )W (2) (x 2 , x 4 )iy (2) (x 3 , x 4 ) 
+^ (2) (x 1 ,x 4 )^ (2) (x 2 ,x 3 )^ (2) (x 2 ,x 4 )^ {2) (x 3 ,x 4 ) (wf^x^+wf^x,)) 
+^ {2) (x 1 ,x 2 )^ (2) (x 1 ,x 3 )^ {2) (x 2 ,x 4 )^ {2) (x 3 ,x 4 ) (y (2) ( Xl ,x 4 ) + ^ {2) (x 2 ,x 3 )) 
+^o (2) (x 2 ,x 3 )^ (2) (x 2 ,x 4 )^ (2) (x 3 ,x 4 ) (w (2) ( Xl ,x 2 ) + ^ (2) ( Xl ,x 3 ) + ^ {2) ( Xl ,x 4 )) 
+W (2) ( Xl ,x 4 )iy (2) (x 2 ,x 4 )^ (2) (x 3 ,x 4 ) (iy (2) ( Xl ,x 2 ) + vy (2) (x 1 ,x 3 )) 

+ (iy (2) (x 2 ,x 4 ) + ^ (2) (x 3 ,x 4 )) (iy (2) (x 1 ,x 2 )iy (2) (x 1 ,x 4 )iy (2) (x 2 ,x 3 )+ 

+l^ (2) (x 2 , x 3 )^ (2) ( Xl , x 4 )W (2) (x 1; x 3 ) + ^ (2) (x 1; x 3 )l^ (2) (x 2 , x 3 )^ (2) ( Xl , x 2 )) 
+iy (2) (x 1 ,x 3 )^ (2) (x 1 ,x 4 )^ (2) (x 1 ,x 2 ) (V (2) (x 2 ,x 3 ) + I^ (2) (x 2 ,x 4 ) + ^ (2) (x 3 ,x 4 )) 
+W (2) ( Xl , x 2 )^ (2) (x 1; x 3 )l^ (2) (x 2 , x 4 )I^ (2) (x 3 , x 4 ) 

+iy (2) (x 1 ,x 2 )^ (2) (x 1 ,x 3 ) (wf } ( Xl ,x 4 ) + iy (2) (x 2 ,x 4 ) + iy (2) (x 3 ,x 4 )) 
+iy (2) ( Xl ,x 4 )vy (2) (x 2 ,x 3 ) (vy (2) ( Xl ,x 3 ) + iy (2) (x 2 ,x 4 ) + iy (2) (x 3 ,x 4 )) 
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+iy (2) (x 2 ,x 4 )^ (2) (x3,x 4 ) (iy (2) ( Xl ,x 2 ) + iy (2) (xi,x 3 ) + iy (2) (x l5 x 4 ; y 
+ (iy (2) ( Xl ,x 2 ) + m/ (2) ( Xi ,x 3 )) (iy (2) (x 2 ,x 3 )iy (2) (x 2 ,x 4 ) + iy (2) (x 2 ,x 3 )iy (2) (x 3 ,x 4 ; 

+ W (2) (X 1 ,X 2 )^ (2) (X 1 ,X 4 )^ (2) (X 3 ,X 4 ) + ^ (2) ( Xl , X3 )W (2) ( Xl , X 4)Wo r 



x 2 ,x 4 ) . 



(65) 



On the other hand, for generic momenta configurations and for models for which onl and &nl 
are sizeable we can expand the exponential entering in the expression (jHHjl of the temperature 
anisotropies up to third order and use the kernel JC(Q in Eq. f)fi()j) . In this way we are able 

re trispectrum of the CMB 



id ]. The parameter g NL is 



to determine the non-linearity parameter which enters into t 
anisotropies according, for example, to the analysis of Refs. (si 
defined through the expansion of the (Bardeen) gravitational potential $ up to third-order 

as 

$ = $L + /NL * ($L) 2 + SNL * ($L) 3 , (66) 

where $l = — 0i, is the linear Gaussian part of $. We find the following expression 

25 25 25 5 

# NL (ki, k 2 , k 3 ) = — (6 NL - 1) + — (a NL - l)-4(k 1; k 2 , k 3 ) + yC(ki, k 2 , k 3 ) - - (a NL - 1) 



+ 



54 



(k 1 -(k 1 + k 2 ))(k 2 -(k 1 + k 2 )) 1 k x k 2 



|ki + k 2 | 4 



3 |k x + k 2 | 2 



cycl. 



(67) 



where ^4(k 4 ,k 2 ,k 3 ) and C(k 1; k 2 ,k 3 ) are defined through Eqs. and (jHSJ- 



V. CONCLUSIONS 



In this paper we showed how to calculate exactly the n-point correlation function of CMB 
anisotropies in the case in which all wavelengths are beyond the horizon at last scatter- 
ing. In this limit the Sachs- Wolfe effect is predominant and its contribution to higher-order 
correlation functions yields the most direct signal of non-Gaussianity in the primordial cos- 
mological seeds. Of particular interest are the bispectrum and the trispectrum which can 
be used to assess the level of non-Gaussianity on cosmological scales. We have calculated 
the non-perturbative expressions for the bispectrum and the trispecturm as predicted within 
single-field models of inflation and in the so-called "squeezed" limit in which some of the 
wavenumbers are much smaller than the others. For other scenarios of generation of the 
cosmological perturbations, we have provided the non-linearity parameters /nl an d #nl 
entering respectively the theoretical predictions of the bispectrum and trispectrum. Our 
results for the bispectrum and the trispectrum represent the essential input in order to 
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obtain the predicted angular bispectrum Bi t i 2 i 3 = (ai imi ai 2m2 ai 3m3 ) and the trispectrum of 
CMB anisotropies (and higher order correlation functions) according, for example, to the 



formalism developed in Refs. 3J, |35|. In these works it is shown how to compute the an- 



gular bispectrum accounting for a non-trivial wavenumber dependence of the non-linearity 
parameter / NL (k 1 ,k 2 ). The angular modulation of the quadratic non-linearity predicted by 
Eq. (|59*|) is currently under investigation 3a ]. adopting the technique of Ref. 3^|, in order 
to look for specific signatures of inflationary non-Gaussianity in the CMB. Notice that for 
a ACDM cosmology a late integrated Sachs- Wolfe effect arising from the explicit time de- 
pendence of the linear gravitational potential during the late accelerated phase would also 
give a contribution on large scales. The formalism developed in this paper can be extended 
to take into account also for this effect. On the other hand on smaller scales there will be 
other effects contributing to CMB non-Gaussianity such as gravitational lensing, Shapiro 
time-delays and Rees-Sciama effects produced at the non-linear level. One should be able 
to compute the angular connected correlation functions induced by these effects by using 
the techinique developed in Ref. 22J and to distinguish them from the large-scale Sachs- 
Wolfe effect provided here thanks to their specific angular dependence. Our predictions for 
the higher-order correlation functions should be compared model by model with the un- 
avoidable contributions from various secondary anisotropies and systematic effects, such as 
astrophysical foregrounds. 



VI. APPENDIX 

In this Appendix we will show in detail how to compute the n-point connected correlation 
functions (J53j) by making use of the generating functional u>(x 1; x n ) and in particular 
how to get the result (|5T|). 
Let us consider the quantity 

5 T = e v+K(ip) - 1 , (A.l) 

where p is a generic Gaussian random field, and K(p) is a generic functional of (p (apart 
from multiplicative coefficients <p will be identified with the comoving curvature perturbation 
( and K will be given by the kernel /C through the iterative procedure). 
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The generating functional for the correlated functions of St is given by 

Z[J\ = [ V^VWl^^^-V , (A.2) 



where J(x) is an arbitrary external source. The functional integral is over all the (p config- 
urations weighted by the Gaussian probability density functional 

e~5 / d y dx ¥ , (y)S(y' x )v( x ) 
V ^ ~ jv[ip}e-^ I ' < ^ A ' 3 ' ) 

which has been properly normalized in such a way that the total probability equals unity, 
fV[<p]V[<p] = l. 

To compute the functional derivatives with respect to J we find it convenient to use an 
additional arbitrary source A(x). We introduce the following generating functional 

Z[J, A] = J p^]p^] e i/^W(e- +,f ^-l) e ^dxA(x) v (x) ; (A 4) 

which reduces to the expression (|A.2|) when A = 0. Functionals of the form in Eq. (|A.4|) 
are common in field theory when computing correlation functions of c omp osite operators (in 
which case J(x) represents a "local coupling"; see for example Ref. |2fij). The correlation 
functions generated by Z[J, A] are given by 

(( e v(yi)+*(*(yO) _ i)...( e v(yn)+*-(v(y»)) _ i)^( Xl )...^(x m )) 

S n+m Z[J,\] 



i 



6J(y 1 )...6J(y n )6X(x 1 )....5X(x m ) 



(J,A)=0 

(A.5) 



Thus we will take only derivatives with respect to J to get the correlation functions of 

e <fi+K(ifi)^ jf we now wr j£ e 

00 K n (ir)\ 

e <P+K( v ) _ 1 = g tL-V£l e <p + ( e v _ 1) ; (A.6) 

n=l n - 

we can rewrite Eq. (|A.4j) as 



Z[J,A]=e J V ' /Z [J,\], (A.7) 

where Z [J, A] corresponds to the generating functional Z[J, A] when the function K{ip) = 

Z [J,X] = Jv[ip}V[ip}e i f d * Ji ^- 1) e i f d * x( * M * ) 

= e Wo[J ' x] . (A.8) 
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In writing Eq. (|A.7|) we have made use of the property i 1 5e l $ dyX ^ v ^ / 5J(y) 
(p(x.)e l f dyX( > y,>ip ( y \ in order to isolate the "interaction term" (see, for example, Ref. \2t\ ) 



Now let us write 



W[J, A] = In Z[J, A] = W [J, A] + In 1 + e 



-Wo 




e 



J dxj(x) £ 



n=l 




) 



(A.9) 



V 



The derivatives with respect to J evaluated for (J, A) = will give the connected correlation 
functions we are looking for, accounting also for the kernel K. We are using the standard 
procedure to evaluate the connected correlation functions, for example, for an interacting 
scalar field (see e.g. Ref. |3])> except that in our case the "interaction term" is related to the 
kernel K and we will make a perturbative expansion around Wo [J, A], since the derivatives of 
W [J, A] with respect to J (evaluated for (J, A) = 0) give the connected correlated functions 
for (e* - 1), see Eq. flU}. 

At this point we have to perform a perturbative expansion. We suppose that the pertur- 
bation parameter is the small r.m.s amplitude of the perturbations, y? rms <C 1 and use the 
fact that the function K((p) is obtained through the iterative procedure described in Section 
III. Let us suppose that 



where the star denotes a convolution operation in configuration space. For the sake of 
simplicity, in the following we will neglect the convolutions, and treat a,b,.. as constant 
coefficients. In Eq. (|A.9J) we will focus on the term ln(l + A) where the definition of A is 
obtained by comparison with Eq. |A~9|) . Thus ln(l + A) = 1 + A - A 2 /2 + A 3 /3 + .... 
Now let us consider how to compute, for example, the three-point correlated function given 
in Eq. (|5ip. At lowest order in our approximation, the only term that in this case is relevant 
is just A and moreover we keep only the first term coming from the expansion of the 
exponential 



Here we have to expand once more, and looking at Eq. (jA.10|) we just need for the bispectrum 
to take K = atp 2 , and from (eiix = 1 + 8/(i5X) + ....) we just pick up the factor 1.. 
Thus we are left with 



K(if) = a*(p 2 + b*ip 3 + .... , 



(A.10) 




(A.ll) 



ln(l + A) ~ A ~ iae 
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i l a J dxj(x) 



' 5W 



+ 



S 2 W 
5A 2 (x) 



(A.12) 



Therefore the bispectrum is given by 

^ {3) (yi,y 2 ,y 3 



S 3 W[J,X] 



+i 



SJ(yi)5J(y 2 )5J(y 3 ) 
(( e ^(yi) _ i)( e ^(y 2 ) _ i)( e *>(ys) _ 1))conn 
5 3 A[J,A] 



(J,A)=0 



(A. 13) 



(J,A)=0 



^(yi)5J(y 2 )5J(y 3 ) 

where we have to compute the last contribution in Eq. ()A.13|) . We need an expression for 
Wq[J, A]. Since the derivatives of Wq with respect to J and A (evaluated in (J, A) = 0) give 
the connected correlation functions for (e^ — 1) and ip, we can write 



oo { n 



n=l 



nl 



W [J, A] = J2~i dxi....dx n z(xi....x n ) J(xi)...J(x n ) 



(A. 14) 



where J(xj) can be either the source J(xj) or A(xj), and correspondingly ■u> n (x 1 ....x n ) are 
the connected correlation functions. For example ((e 1 ^* 1 ) — I)<y2(x2)v?(x3)) for the choice 
J(x!)A(x 2 )A(x3), and one has to consider different combinations. 



Using the usual operations for functional derivatives (see, for example, 27]) we find 



5 3 A[J,X] 



SJ(yi)5J(y 2 )5J(y 3 



= a J2[w2 (y P1 , y P2 w 2 (y P3 , y P2 ) + -z 4 (y Pl , y P2 , y P3 , y P3 

(J,A)=0 V Z 



(A.15) 

where the sum is over the permutations pi,p 2 ,P3 of indices (1,2,3) and we have used the 
following notations 



tf5 2 (x,y) = ((e^ - lMy)> conn . 
«; 4 ( Xl ,x 2 ,x,x) ee ((e^)-l)(e^-lV(x)^(x)) c 



(A.16) 



Notice that it is possible to compute the connected correlation functions appearing in 
Eq. |A~Tfijl and ifOTjl in a similar manner to what one does for e^. 

The result in Eq. (IA.15J) is a sum of two terms that correspond to the two pieces in 
Eq. ()A.12j) . One has two evaluate <p at two equal points because we have taken the derivative 
twice with respect to A. Finally the product of the two-point connected correlation functions 
and the presence of the fourth-order connected correlation function are due to the fact that 
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we have the product of two first derivatives w.r.t A and a second order derivative w.r.t A, 
respectively. Taking then the three derivatives w.r.t. J involves (e v — 1) in the correlations. 
Now it is easy to generalize the result (|A.15|) when a is not a constant coefficient but we 
have a kernel in configuration space such that 

K(ip) = J dx 1 dx 2 K(-x — X!,x — x 3 )</?(x 1 )<^(x 2 ) + .... (A. 17) 

As one can guess from Eq. (|A.15|) one has 



t- 3 



5 3 A[J,X] 



(J,A)=0 



SJ(yi)5J(y 2 )5J(y 3 ) 
E p J rfx!rfx2i^(y P2 - xi,yp 2 - x 2 )w 2 (yp 1 ,x 1 )w 2 (yp3,X2) 

+ \ll J dxicZx 2 K(y P3 - xi, y P3 - x 2 )w 4 (y Pl , y P2 , x l5 x 2 ) 



(A.18) 



In fact we have explicitly verified that this is the correct generalization of Eq. (|A.15|) . 

Thus from Eq. (jA.9|) and Eq. (jA.18|) the connected three-point correlation function for 

( e <p+K(<p) _ readg 

S 3 Z[J,\] 



W ( V(y 1 ,y 2 ,y 3 )=i- 3 



(J,A)=0 



(A.19) 



SJ(yi)5J(y 2 )5J(y 3 ) 
( (e ^(yi) _ i)( e ^(y 2 ) _ i)( e *>(y 3 ) _ i)) connected 

+ Yj dx 1 d^ 2 K(y P2 -5c 1 ,y P2 -x 2 )w 2 (y pi ,^ 1 )iv2(yp 3 ,^ 2 ) 
p 

+ \Y J dxidx 2 fT(yp3 -xi,y P3 - x 2 )w 4 (yp 1 , y P2 , xi, x 2 ) , 

where one has to use the definitions in Eqs. (jA.16|) and (jA.17|) . 
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